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Abstract. We investigate the dynamic impact of heterogeneous environments on 
superdiffusive random walks known as Levy flights. We devote particular attention to 
the relative weight of source and target locations on the rates for spatial displacements 
of the random walk. Unlike ordinary random walks which are slowed down for all 
values of the relative weight of source and target, non-local superdiffusive processes 
show distinct regimes of attenuation and acceleration for increased source and target 
weight, respectively. Consequently, spatial inhomogeneities can facilitate the spread of 
superdiffusive processes, in contrast to common belief that external disorder generally 
slows down stochastic processes. Our results are based on a novel type of fractional 
Fokker-Planck equation which we investigate numerically and by perturbation theory 
for weak disorder. 
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A particle performing ordinary diffusion is typically characterized by spatiotemporal 
scaling relation ~ t 1 / 2 . An increasing number of physical and biological 

systems are in conflict with this relation and exhibit anomalous diffusion. Whenever 
~ t 1 ^ with an exponent < /i < 2 a system is said to exhibit superdiffusive 
behavior. Superdiffusion was discovered in a wide range of systems, for instance chaotic 
dynamical systems [I], particles in turbulent flows [2j, saccadic eye movements [3], 
trajectories of foraging animals [4], and very recently in the geographic dispersal of 
bank notes [5]. 

Theoretically, superdiffusion is often accounted for by scale free random walks 
known as Levy flights for which successive spatial displacements Ax are drawn from 
a probability density function (pdf) with an algebraic tail, i.e. p(Ax) ~ |Aa;| _d_ ^, 
where d is the spatial dimension and the Levy exponent /x < 2, such that the variance of 
Ax is divergent [6]. Embedded in the context of continuous time random walks (CTRW) 
one can show that the dynamics of the pdf of finding a superdiffusive particle at position 
x at time t is governed by a fractional generalization of the diffusion equation, 

d t p(x,t) = D,A^ 2 p(x,t), (f) 

in which the ordinary Laplacian A is replaced by the fractional operator A^/ 2 , a non- 
local integral operator, the action of which is equivalent to a multiplication by — \k\^ in 
Fourier space. The constant is the generalized diffusion coefficient. 
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A process governed by Eq. (|TJ) is spatially homogeneous and isotropic, as the 
probability rate w(y\x) of a displacement from x to y depends only on distance, i.e. 
w(y\x) oc \x — y\-( d+ ^. 

Numerous random processes, however, occur on spatially disordered substrates 
or evolve in the presence of quenched spatial inhomogeneities. Depending on the 
underlying physical or biological system, one obtains various generalizations of (pQ) 
which incorporate the spatial structure. In the generalized Langevin approach [7], [8], [9 J an 
additional force term — VF on the right hand side of (pQ) accounts for an external position 
dependent force field F(x). In topologically superdiffusive systems, such as intersegment 
transfer of gene regulatory enzymes on DNA strands [TQl CEH CE2] , the transition rate is 
modified by a Boltzmann factor, i.e. w(y\x) oc \x — y x exp —f3[V(y) — V(x)]/2, 

where V(x) is a position dependent potential and (3 an inverse temperature. In 
subordinated superdiffusive processes, an ordinary diffusion process subjected to an 
external force is sampled at highly variable operational time intervals [13]. All three 
systems exhibit very different response properties to the imposed spatial structure [14], 
but converge to the same Fokker-Planck equation in the limit of ordinary diffusion. 

Here we investigate the dynamic impact of heterogeneous environments and devote 
particular attention to the relative impact of source and target locations on the rates 
w(x\y) for spatial displacements of the random walk (see FigH]). We define the spatial 
inhomogeneity in terms of the attractivity or salience s(x) > of a location x. For large 
and small values of s(x), the location x is attractive and unattractive, respectively. We 
assume that in equilibrium a walker's stationary probability p* (a?) of being at a location 
x is proportional to the salience at a?, i.e. p*(x) oc s(x). In that respect, the salience field 
can be defined operationally as the likelyhood of finding a walker at site x. Furthermore 
we assume that a transition from y to x is more likely to occur when the salience is 
large at the target location x, and less likely to occur when the salience is large at the 
source location y. This leads to a transition rate 

w(x\y) = -s c (x)f(\x-y\)s c - 1 (y), (2) 

T 

where r is a time constant and the sandwiched term f(\x — y\) oc \x — with 
< /i < 2 accounts for Levy flight jump lengths. Inserting j2]) into the master equation 

dtp(x,t) = J dyw(x\y)p(y,t) - w(y\x)p(y,t) (3) 

one can see that detailed balance is fullfilled and that the stationary solution if it exists 
is proportional to the salience. The central parameter in our analysis is the weight 
parameter < c < 1, which quantifies the relative impact of source and target salience 
on the dynamics. When c = 1 a transition y — > x only depends on the salience at the 
target site and is independent of the salience at the source, when c = the salience of the 
target site has no influence on the transition and the rate is decreased with increasing 
salience at the source. The intermediate case c = 1/2 is equivalent to topological 
superdiffusion approach with a salience given by a Boltzmann factor s(x) = e~^ v< ^ . 
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Figure 1. Random walk processes in inhomogeneous salience fields s(x) in two (a) 
and one (b) dimensions. Source and target locations of a random jump are denoted 
by y and x, respectively. 



In combination with the rate ((2j) the master ((3l) is equivalent to the fractional 
Fokker-Planck equation 



where p = p(cc,t), s = s(x) and D is a generalized diffusion coefficient. The fractional 
Laplacian is defined by 



with = 2> M iT- n/2 T((fi + n)/2)/r(-///2). In Fourier space A^ 2 corresponds to 
multiplication by -|fe|^: ^{A^fjik) = -\ k Y L F{f}( k )- Note that the fractional 
Fokker-Planck (111) is equivalent to a number of known stochastic processes for specific 
choices of the parameters c and \i. For instance when s(k) = const ([4]) is equivalent to 
free superdiffusion, i.e. (DD), when c = 1/2 and n = 2, (pH) reads d t p = —VFp + Ap, 
with F(x) = —V'(x) and V(x) = — logs(cE), i.e. diffusion in an external force 
field. When c = 0, (pj reads d t p = DA^ 2 p/s, which is equivalent to generalized 
multiplicative Langevin dynamics for the process X(t), i.e. dX = D(X)dL fl with 
D(X) = v^2(exp PV (X) / '2) / D and L^(t) is a homogeneous Levy stable process. 

In the following we investigate the relaxation properties of one dimensional 
processes governed by the fractional Fokker-Planck equation ([4]). It is convenient to 
make a transformation of variables ip = s 1 ^ 2 p and rewrite the dynamics as a generalized 
Schri^inger equation d t ip = T~Cip with a Hamiltonian 



possessing identical spectral properties. 

We model inhomogeneities as salience fields of the type s = exp(— ev), where 
v — v(x) is a random potential with zero mean and unit variance. The parameter 
e > quantifies the strength of the inhomogeneity. Typically, inhomogeneities posses a 
spatial correlation length. For v we chose random phase potentials defined by 



d t p = D s c A» /2 s {c - 1] p -Dp s^A^s", 



(4) 




Hip = s 



72^/2^-1/2^ _ ^ s c- lAM /2 s c 



(5) 




(6) 
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Here the phase 8(k) is uniformly distributed on the interval (0, 2ir] and the amplitude 
is given by the correlation spectrum S(k), 4>{k)(f){k') = 2nS(k)5(k — k'), which we 
chose to be Gaussian with a correlation length £: S(k) = 2£exp(— A; 2 £ 2 /7r). For small 
e we can expand ([5]) up to the second order and obtain a Hamiltonian of the form 
Ji = A^ 2 + U which we can treat perturbation theoretically to obtain a spectrum of 
the form E{k) ~ —D^ c (k;e)\k\^. Here D^ c (k;e) describes the relaxation properties 
on the corresponding scale length A = k~ x . In the limit of vanishing potential the 
generalized diffusion coefficient is identical to the diffusion coefficient D of the free 
system. Typically, a random spatial inhomogeneity slows down a random process, and 
one generally expects D^ c {k;e) to be smaller than D^ c (k;e = 0). Up to the second 
order in e we obtain 

D tl , c (k;s)/D = l-e 2 G tl , c (k), (7) 

where G^ c (k) = ±f dqS(q)g^ c (k/q) and 

g,, c (z) = - 2(c - 1/2) 2 + ^r{2c(c - 1) - (c - 1/2) 2 (|1 - z\» + |1 + z\») 

+ [(c - l/2)(z» + |1 + z\») - c] 2 /(|l + z\" - z") 

+ [(c - l/2)(^ + |1 - z\») - c] 2 /(|l - z\» - (8) 

The limit k — > yields the long scale asymptotics for the process for which we 
obtain 

G m _ < f 1 /2-c 2 ,0< / ,<2, 

^ c(0) -\ 1/2 + 2^,^ = 2, (9) 

This result indicates that with increasing influence of the target salience (increasing 
values of c) superdiffusive processes (// < 2) exhibit the opposite behavior as ordinary 
diffusion processes (// = 2). As c is increased G^ c (0) increases as well for ordinary 
diffusion, which means that these processes are attenuated more strongly. Quite contrary 
to superdiffusive processes, for which a more pronounced target influence decreases 
G^ c (0). This implies that slowing down of superdiffusive processes becomes weaker 
as the target weight in the rates is increased. Note also that for weak potentials the 
magnitude of this acceleration is independent of the Levy exponent /i, see FigfJl Only 
when c = 0, and target salience has no impact on the transition rate, all processes are 
slowed down by the same amount, i.e. G^^O) = 1/2. 

Note that, as c is increased for \i < 2, the function G^ c (0) even becomes negative 
for c > c cr it = 1/a/2. This implies that these processes are no longer slowed down by the 
spatial inhomogeneity but rather accelerated, as a negative value for G MiC (0) implies a 
diffusion coefficient larger than the one for free superdiffusion, i.e. D^ c {k — > 0; e) > D. 
Consequently, the common notion that random processes are typically slowed down by 
spatial inhomogeneities is not valid when superdiffusive processes are involved. 

For large e we computed the generalized diffusion coefficient D^ie) numerically for 
three periodic and one random phase potential: a cosine potential v(x) = \/2cos(x/A), 
a potential with localized potential minima and one with localized potential peaks, for 
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Figure 2. The impact of source and target location on the asymtotics of random walk 
processes. The generalized diffusion coefficient depends quadratically on the magnitude 
eof the inhomogeneity: D^ tC = 1 — G MiC (0)e 2 . The pre-factor G M)C (0) quantifies the 
impact of the spatial inhomogeneity, i.e. when G M , C (0) > the process is slowed down, 
when G AtiC (0) < the process is accelerated. The two lines depict G AtjC (0) as a function 
of the weight parameter c for ordinary diffusion (fx = 2, upper line) and Levy flights 
(/i < 2, lower line). As c is increased (increasing weight of the target location), ordinary 
diffusion is slowed down more strongly contrary to Levy flights for which attenuation 
decreases until a critical value of c cr it = l/v2 is reached (denoted by small circle). 
Beyond this point Levy flights are accelerated by the external inhomogeneity. 

which v(x) = ±acos 7 (x/A) + b (a, b > 0, 7 = 32), and a potential defined by (E|) with a 
gaussian spectrum. The results are depicted in FigJH 

For c = (i.e. full weight of source location) the processes in a given environment 
behave similarly, independent of the exponent \i (Figj3](a)), ordinary diffusion and all 
superdiffusive processes exhibit the same quadratic decrease of D^ c with e in a fixed 
potential. For c = 1/2 (Figj3](b)) only the superdiffusive processes exhibit an identical 
response to a given potential, the response of the ordinary diffusion process differs. When 
c is increased beyond the critical value the difference between ordinary diffusion and 
Levy flights becomes maximal and changes qualitatively (FigHj). In this regime ordinary 
diffusion processes are still slowed down by the spatial inhomogeneity as opposed to the 
Levy flights which exhibit an increase in the diffusion coefficient with increasing potential 
strength and are thus accelerated. For small values of e the numerics agree well with our 
results obtained by the perturbation theory, above (i.e. (jnj)). Despite the fact, that for 
all values of /x < 2, the curves collapse on one single function, Levy flights are sensitive 
to the potential shape. The relatively pure performance of the random phase potential 
(D on the Figj3)) is due to the finite size effects, which does not constitute an obstacle for 
the periodical potentials (A,B,C). Therefore we have omitted the results for the random 
phase potential in case c = 1 (Fig. 0J). 

The above considerations were restricted solely to the infinite systems. A key 
question is how these processes behave in finite systems and to what extent finite size 
effects play a role. Therefore we investigate the relaxation properties in finite systems 
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Figure 3. Generalized diffusion coefficient Z3 M:C as a function of the magnitude e of 
the salience field six) — exp[— ev(x)]. The function Z? MjC is shown for four types of 
potentials, a cosine potential (A), potentials with localized maxima and minima, (B 
and C, respectively) and a random phase potential (D). Potentials are displayed on 
the top panel. Results for two values of the weight parameter c are presented on the 
panel (a) c = and (b) c = 0.5 correspondingly. Dashed line and dash-dotted lines are 
analytical results for [i — 2 and < /i < 2 we obtained from perturbation theory ((9j). 
Each type of symbol corresponds to the one Levy exponent: O - fx = 0.5, <)-/-* = 1-0, 
A - fi= 1.5, □ - fi= 2.0. 




Figure 4. Generalized diffusion coefficient D^ tC as a function of the magnitude e of 
the different salience fields for the case c—1. The meaning of the symbols see caption 
to FigO In contrast to diffusive processes, which are attenuated, the superdiffusive 
ones are enhanced by inhomogeneity. 
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Figure 5. The impact of the heterogeneity on the relaxation time of the processes 
St ~ g^ c (X/L) with different power-law exponents fj,, in dependence on the size of the 
system L for the cosine potential with period A. Gray lines denote the large scale limit, 
according to ([9]) . 

of length L. To this end we consider the quantity 

St = r/r frcc - 1, 

where r is the relaxation time of the process and rf ree is the relaxation time of 
the same process without spatial inhomogeneity. For the cosine potential we obtain 
5t « e 2 g^ c (X/L), where g^ c is defined in (JHJ). See Figj5] for three values of c = 0, 0.5 
and 1. Small values of 5t correspond to a small effect of the inhomogeneity. For c = 
(Fig.[5](a)), the ordinary diffusion process exhibits the smallest value of 5t, which implies 
that in situations in which the source salience is important, diffusion relaxes fastest. On 
the contrary, for c = 1 (Fig. [5](c)) strongly superdiffusive processes (i.e. \i — > 0) exhibit 
a small St. Only when both, source and target possess an equal impact on the jump 
rates (i.e. c = 0.5 j5](b)), 6t exhibits a minimum for intermediate values of the Levy 
exponent /i [15]. 

We considered the consequences of the relative weight of source and target locations 
in one-dimensional random walk processes, evolving in inhomogeneous environment. 
Our analysis revealed essential differences between superdiffusive Levy flights and 
ordinary random walks when they occur in regular and random spatial inhomogeneities. 
Unlike ordinary random walks, Levy flights can be accelerated when the influence 
of the target salience is sufficiently large, which may shed a new light on optimal 
search strategies in heterogeneous landscapes, and dispersal phenomena in population 
dynamical systems and various physical and biological contexts. 
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